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Abstract. We have theoretically investigated the focusing of a launched cloud of 
cold atoms. Time-dependent spatially-varying magnetic fields are used to impart 
impulses leading to a three-dimensional focus of the launched cloud. We discuss 
possible coil arrangements for a new focusing regime: isotropic 3D focusing of 
atoms with a single magnetic lens. We investigate focusing aberrations and find 
that, for typical experimental parameters, the widely used assumption of a purely 
harmonic lens is often inaccurate. The baseball lens offers the best possibility for 
isotropically focusing a cloud of weak-field-seeking atoms in 3D. 



1. Introduction 

The development of laser-cooling techniques pQ has facilitated the preparation of 
samples of atoms at microKelvin temperatures 0- It is therefore now possible to 
modify drastically the centre-of-mass motion of atoms, in direct contrast with the small 
angular deflection of fast beams studied prior to the development of laser cooling 3 . 
The (very strong) optical scattering force utilised to cool atoms can also be used 
to manipulate their position, as can the (relatively weak) conservative optical dipole 
force. Since the dipole force is coherent it represents a better way of manipulating 
cooled atoms than the scattering force (for which the random nature of absorption 
and spontaneous emission leads to Brownian heating). An alternative method of 
manipulating paramagnetic cold atoms is to use the Stern-Gerlach force 0]. 

One of the goals in the field of atom optics [5] is to realise atom-optical elements 
that are analogues of conventional optical devices, such as mirrors, lenses and beam- 
splitters. An atom mirror reverses the component of velocity perpendicular to the 
surface and maintains the component parallel to the surface. An atom lens can 
modify both the transverse velocity component and the longitudinal component. To 
date, the Stern-Gerlach force has been used to realise flat atomic mirrors curved 
atomic mirrors [3 , and pulsed mirrors for both cold (thermal) 8 and Bose condensed 
atoms 9 . It has also been demonstrated that the surface of a magnetic mirror can 
be adapted in real time with corrugations that can be manipulated in times shorter 
than the atom- mirror interaction time |1U|. 

There are many reasons for studying focusing of launched cold atoms: atom 
lithography JI] ; transferring cold atoms from a MOT to a remote vacuum chamber of 
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lower background pressure |12j ; loading miniature magnetic guides |13| , atom chips | 14| 
and storage rings . In comparison to an unfocused cloud, the density of the cloud 
can be increased by many orders of magnitude after magnetic focusing. 

The first demonstration of 3D focusing using pulsed magnetic lenses was 
conducted by Cornell et al. ^B]. However, their work did not address the optimum 
strategy for achieving a compact focused cloud, nor did they discuss the limiting 
features for the quality of their atom-optical elements. The group of Gorceix has 
performed experiments demonstrating the longitudinal Stern-Gerlach effect with an 
atomic cloud using pulsed magnetic forces |17|. and an experimental and theoretical 
study of cold atom imaging by means of magnetic forces |18j. To date, the experimental 
and theoretical studies of pulsed magnetic focusing have been analysed under the 
assumption that the magnetic lens potential is harmonic - this work addresses the 
validity of this approximation, and the effects of aberrations. We note also that there 
has been no theoretical or experimental work on full 3D focusing of cold atoms using 
a single magnetic pulse, although this is a very useful focusing strategy. 

The scope of this paper is to investigate theoretically and numerically the 
limiting factors to the quality and size of the final image obtained in pulsed magnetic 
focusing experiments; to identify the sources of aberration; and to discuss schemes for 
minimising their deleterious effect. In this paper we restrict our attention to focusing 
strategies using a single magnetic impulse. A second paper describing more than one 
impulse (alternate gradient focusing) is in preparation. Whilst this work shall 
concentrate on the analysis for achieving a compact cloud in space, it is also possible 
to use pulsed magnetic fields to reduce the momentum spread of an expanding cloud 
with appropriate magnetic impulses. This can be viewed as an implementation of 5- 
kick cooling, which has been demonstrated with atoms [20], ions and Bose Einstein 
condensates (BEC) ,9> Some of the techniques described here are also successfully 
used for the deceleration of polar molecules using time- varying electric fields |22- 
Atom-optical elements realised with light forces J] EH] are beyond the scope of this 
paper. 

The remainder of the paper is organised as follows: Section 2 outlines the theory 
of how to achieve the desired magnetic fields; Section 3 contains an analysis of 
magnetic imaging and minimising the final cloud size; Section 4 describes and contrasts 
the performance of different magnetic lenses; Section 5 contains a discussion and 
concluding remarks. 

2. Magnetic lens theory 

2.1. The Stern-Gerlach force 

An atom's magnetic dipole interaction energy is U = mpgp/iBB for a field of 
magnitude B, where mp is the magnetic quantum number, gp is the Lande g-factor 
and [1b the Bohr magneton. We assume that field zeros are avoided and that the 
magnetic moment adiabatically follows the field. Depending on an atom's hyperfine 
quantum state, the Stern-Gerlach force F = —VU cx VB can be used to attract it 
towards weak (if mpgp > 0) or strong (if mpgp < 0) magnetic fields. The choice 
of whether atoms in weak or strong-field seeking states are launched depends on the 
particular application. For some applications, e.g. loading a remote dipole trap, or 
a secondary magneto-optical trap, it does not matter which atomic state is used. In 
this work we assume that the ensemble (of alkali metal atoms) has been optically 
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pumped into the stretched state \F = I + 1/2, rap — F) (a weak-field seeking state 
with mpgp = 1). 

A purely harmonic magnetic field magnitude will result in an aberration-free lens. 
In this paper we will consider current-carrying wires assembled to give two different 
kinds of second order magnetic field magnitude: 

B 1B (x, y lZ ) = B Q + ^ (~x 2 /2 - y 2 /2 + (z - z c f) , (1) 

B 3D {x, y, z) = B + B 1 {z - z c ) + -y (x 2 + y 2 + (z - z c ) 2 ) . (2) 

Bo, B\ and B2 are the bias field, the axial gradient and the field curvature, respectively. 
In both cases these parameters are chosen to prevent the field magnitude exhibiting 
zeros in the region of interest, thus avoiding Majorana spin-flip transitions [23]. The 
point {0, 0, z c } defines the centre of the lens. A lens of the form of -Bid can be used 
either to focus axially or radially; a lens of the form of B3D is used to focus isotropically 
in 3D. The accelerations associated with these fields are harmonic about {0, 0, z c } : 

ai D = -tu 2 {-x/2,-y/2,(z- z c )}, (3) 

a 3 D = -u 2 {x,y, (z - z c )} + a , (4) 

where uj 2 — ii^B2/m is a useful measure of the power of the lens, ao = {0, 0, /iB-Bi/rre} 
is a constant acceleration arising from axial magnetic gradients, and m is the atomic 
mass. Note that for the ID case (equation (J3J) the axial curvature is twice the 
magnitude of and opposite in sign to the radial curvature, ui 2 — —2io 2 . 

The acceleration of equation J3) results in a lens which is axially converging 
(diverging) and radially diverging (converging) if B2 is positive (negative). In order to 
achieve a 3D focus with such lenses, an axially converging lens pulse must be followed 
by an appropriately timed axially diverging lens (or vice versa). This alternate- 
gradient focusing strategy is the subject of another paper |19| . In contrast, this paper 
details novel lenses which yield accelerations of the form in equation |@} . Isotropic 3D 
focusing can be achieved with a single lens pulse for atoms in weak-field-seeking states, 
as Earnshaw's theorem (which states that a maximum of magnetic field magnitude in 
free space is not allowed [25]) ensures B2 > in this case. Note that the harmonic 
accelerations in both equation © and lead to three separable one-dimensional 
simple harmonic oscillator equations for the atomic motion. 

2. 2. Magnetic fields from current bars and circular coils 

This work considers straight current bars and circular coils for the formation of 
lenses. The Biot-Savart law yields magnetic fields that are analytic for both finitc- 
and infinite-length current bars. For circular coils the field can be expressed in 
terms of elliptic integrals [23 • A discussion of the form of the contours of magnetic 
field magnitude for various magnetic trapping configurations has been provided by 
Bcrgcman et al. [2Z|. The fields are constrained by Maxwell's equations, which, in 
conjunction with symmetry arguments, allow the spatial dependence of the fields to 
be parameterised with a small number of terms. In particular, for a cylindrically 
symmetric magnetic coil configuration, the fourth-order on-axis ID magnetic field, 
B z (r = 0, z) = Y^,t=o Bi z l /il, gives the complete fourth-order 3D magnetic field: 

B(r,z)={B r ,B z } 
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--r--rz,B + Bl z + -[z 2 

{r 3 rz 2 z 3 r 2 zl fr 3 z rz" 2f , . 

with third-order magnitude: 



B(r, Z )=B + B lZ + -B^ + -^—--)r 
2 f B 3 B\ B 2 Bi 3 



If, in addition, the coil system is axially-symmetric the magnitude to fourth order is: 
B(r,.) = Bo + - rV2) + + + + §l) . (7 ) 

Let us consider two coils of N turns with radius a, separation s, carrying currents 
Ii and 12- It is convenient to partition the currents in each coil as a current In with 
the same sense and a current /ah in opposite senses, i.e. 2ijj = I1+I2, 21 An = I\ — l2- 
If one defines r\ = /iqNI/2, and uses the scaled separation S = s/a, the axial magnetic 
field is thus: 

BCfl ^ = f (?7H + ??ah) (r; H - ??ah) \ fR v 

{ , ) \a(l + (z/a-S/2) 2 ) 3 / 2 a{l + {z/a + S/2) 2 ) 3 / 2 J 1 {> 

yielding the axial Taylor expansion terms: 

2 m r 3 Vah S 6 m (S 2 -l) 



a(I + S 2 /4) 5 a 2 {l + S 2 /A)* ' a 3 (l + S 2 /4)^ 

15^^ (5 2 - 3) 45r ?H (S 4 - 6S 2 + 2) 

±»3 — — ) -04 — n — • l y i 

a 4 (l + 5' 2 /4) 5 a 5 (l + S' 2 /4)- 

Thus axially symmetric lens systems (equation (jTJ): 77^^ = 0, S dd = 0), consisting of 
one or two circular coils have radial and axial strengths: 

2_ „ 2 _ 3^^/(1 + sign(S))(S 2 -l) 
2ma 3 (l + ,S 2 /4) 7/2 

The sign function switches between a single and a double coil. The expressions of 
equation 10 will be utilised in the next section where various configurations of coils 
and bars for realising magnetic lenses are considered. 



2. 3. Configurations for realising magnetic lenses 

There are six distinct coil and current bar configurations used in this paper for 
single-pulse focusing. Figure displays these six strategies, whilst figure [5] shows the 
associated magnetic field information. Strategies I-III lead to the axial/radial focusing 
of equation Q), whereas strategies IV- VI yield isotropic 3D focusing (equation Q). 
Strategies I-V deal with cylindrically symmetric coil arrangements. 
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Figure 1. The six different lens strategies detailed in the text. Strategy I utilises 
the centre of a single coil; Strategies II and III use the geometric centre of a pair 
of coaxial coils (carrying equal currents in the same sense) with separations of 
S"=0.58 and S=2.63 coil radii, respectively; Strategy IV uses the geometric centre 
of a pair of coaxial coils with unequal currents and a relative separation S = \/3; 
Strategy V uses a single coil axially offset to z/a = ± ^2/7 ; Strategy VI uses the 
geometric centre of a Baseball coil with dimensions w = I = 2a combined with a 
coaxial coil pair with 5 = 1. 



2.3.1. Strategy I: a single coil The magnetic field magnitude at the centre of a 
single coil of radius a with NI\ current turns is characterised by the coefficients in 
equation 10 with S — and using I2 — when evaluating r/n. The radial curvature 
is positive, thus, a single coil can be used to radially focus weak-field seeking atoms. 
However this lens has an axial curvature that is negative and twice the magnitude. 
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Distance from Centre of Lens / R, X or Z 

Figure 2. Here we show the respective axial and radial variation of the magnetic 
field magnitude associated with the six lens strategies I- VI of Figure HI R. X 
and Z are the scaled radial, x and axial displacements from the centre of the 
lens, respectively. In Images IV and V the axial gradients in the magnetic field 
magnitude alter the centre-of-mass velocity of an atom cloud but do not lead to 
lensing. Parabolic approximations to the field magnitudes are also given. 



Figure n Strategy I shows the geometry of this lens. 

2. 3. 2. Strategies II, III and IV: a pair of axially- displaced coils For a pair of separated 
coaxial coils, where both coils have equal current with the same sense, an axially 
diverging (converging) lens is realised for separations less (greater) than the coil radius 
(equation I|1U|)). The lens curvature is zero when the separation equals the coil radius 
(the Helmholtz condition for achieving uniform fields). The axial curvature is a factor 
of -2 greater than the radial curvature. 

Since strategies II and III are axially-symmetric configurations there are no third 
order terms. The third term in equation JJJ) is zero when B4 = 0, which occurs 

if S = v / 3~±77 = 0.595 or 2.38. The fourth term in equation Q is zero when 
B 4 B = Bl, which occurs if S = ^(13±VT^) = 0.565 or 2.89. The harmonicity 
of a radial-focusing lens (S < 1 Strategy II) is thus optimised if the relative coil 
separation is S = 0.58, whereas the harmonicity of an axial-focusing lens (S > 1 
Strategy III) is thus optimised if the relative coil separation is S = 2.63. Figure Q II 
and III show the geometry of these optimised lenses. 

Figure^Strategy IV shows the geometry of an isotropic 3D lens formed from two 
coils carrying different currents. From Equation © we see that an isotropic 3D lens 
is formed when B\ 2 = 6B2B0, and this can be re-expressed using equation @ as: 

77ah - ±W8(1-1/S 2 )- (11) 

From equations 1)619(1 we see that some of the third-order terms in the magnetic field 
magnitude are removed by setting S = \/3 and thus B3 = 0. This means that 
the coils carry currents of I\ and I2 = — 0.396/i respectively. Note that the only 
way to remove all third-order terms is to have S — V3 and set B\ 2 — 2B2B0 (i.e. 
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'/ah = ±?7h\/8/3(1 — 1/S 2 ) ) which corresponds to a purely axial lens. Purely radial 
lenses are achieved when Bi = 0, (i.e. S — 1) for non-zero tjah and rjn. 

2.3.3. Strategy V: an axially offset single coil At the centre of a single coil the radial 
and axial curvatures have opposite signs (Strategy I), however, the radial and axial 
curvatures have different z— dependence, (Figure ^ Strategy V). At the two axial 
locations z = ±-y/2/7a the curvatures are equal in magnitude and both positive. 
Therefore an impulse applied to a cloud whose centre-of-mass is at either of these 
positions will lead to isotropic 3D focusing. 

2.3.4- Strategy VI: loffe-Pritchard configuration Ioffe-Pritchard (IP) traps are used 
extensively for atom trapping and are similar to the Ioffe configuration utilised in 
plasma confinement |28|. One of the simplest forms of this trap is with a pair of coaxial 
coils to generate both a bias field, Bq, and an axial curvature Bi. Four current bars 
run parallel to the z-axis, with each bar running through the corner of a square (with 
side length w — Wa) centred on the axis. The magnitude of the current in each bar 
is equal, but neighbouring bars have the opposite current sense. 

This configuration is not strictly radially symmetric, although it is an excellent 
approximation for small displacements from the axis. The four current bars give rise 
to a transverse magnetic field of the form {B x ,B y } = B[{x, —y} close to the axis. 
Here, B[ = A^iqI' / (ttw 2 ) , where I' is the magnitude of the current in a bar, and the 
prime distinguishes it from the axial gradient defined previously. The shape of the 
magnetic field magnitude contours depend on the geometry of the trap and the ratio 
of the coil current to the bar current. Such traps are extensively used in Bose-Einstein 
condensation experiments, where an additional circular coil pair is used primarily to 
lower the bias field B n , yielding an anisotropic harmonic-oscillator potential, with tight 
radial confinement B^ r = B[ 2 / Bq — B2/2, and weaker axial confinement B^ z = B2- 
For the purposes of magnetic focusing an isotropic IP trap is required. This is achieved 
if the following relationship holds among the lens parameters: 

ir= 3 ^V(^w" (12) 

In BEC experiments variants of the IP trap described above are used, since 
infinitely-long current bars are not realisable. An elegant winding pattern is the 
baseball geometry. A baseball coil has an axial curvature in addition to a radial 
gradient, and these quantities can no longer be independently varied. FigureQ]Strategy 
VI shows the cuboidal baseball geometry, where the sides have lengths w = Wa, 
w = Wa and I — La, (I is along the z-axis). The bars carry a current It is 
impossible to realise a 3D isotropic lens solely by adjusting the aspect ratio of the 
coils w/l, and for this reason we add a coaxial pair of coils carrying equal currents in 
the same sense (if tjah 7^ then the equality of the x and y curvatures is broken). 
The magnetic field magnitude has a third-order Taylor expansion: 

B = B '» + [Ik - B -r) <* 2 + * 2) + ¥«" + ( B » - w) «" 2 (13) 

where dashed terms indicate contributions solely from the baseball coil, and the 
subscript f is used when one must add together the Bq, B2 Taylor contributions from 
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the circular coils (equation @) and B' Q , B' 2 from the following baseball coil Taylor 
terms: 

^L\V 2 , ^IL (L 5 + 3L 3 W 2 + 4LW 4 ) 

^0 - , TO , TTr n, /T n , „ =?> B l - 



a(L 2 + W 2 )y/L 2 + 2W 2 a 2 W 2 (L 2 + W 2 Y^L 2 + 2W 2 

, i^Il (6L 6 W 2 + 18L 4 W 4 + 11L 2 W 6 - 5W 8 ) 

a 3 (L 2 + W 2 f(L 2 + 2W 2 y 

(-5L 7 W 2 - 10L 5 W 4 + 11L 3 W 6 + 24LW 8 ) 
B' = —5 — i , '-. 14 

a 4 (L 2 + W 2 ) 4 (L 2 + 2W 2 )t 

The simplest isotropic lens to calculate is when we use circular Helmholtz coils - 
i.e. S = 1, B2 = and the only contribution of the circular coils in equation (|13|l is 
the axial constant field Bo, which we tune to enable lens isotropy. The finite- length 
current bars lead to a non-trivial relationship among the parameters for realising an 
isotropic lens by adjusting the circular coil bias-field and curvature. For the purpose 
of our calculations we use the parameters W = W = L = 2 for a cubic baseball, with 
S = 1 and a relative circular coil current of ?jf = 0.154. 



2-4- The harmonicity of the magnetic lenses 

The previous subsection was a description of the way to achieve a radially converging 
lens with a single coil or a pair of coils; an axially converging lens with a pair of 
displaced coils; and three geometries for achieving an isotropic 3D lens. Aberrations 
are caused by departures of the actual potential from the ideal harmonic potential. 
To measure the degree of departure from harmonicity of a lens, as a function of the 
distance from the lens' centre, we use: 

l a F - a H | , s 

|a H - a | 

where an is the harmonic fit, (equation or Q), to the full Biot-Savart lens 
acceleration ap. The departure from harmonicity is a function of position, and the 
cylindrical co-ordinates R and Z are used to plot e(i?, Z) for different coil systems. 
Note the scaled co-ordinates R and Z are normalised to measure length in units of the 
coil radius, R = r/a and Z — z/a. Figure JHJl shows the departure from harmonicity 
for the six different focusing strategies. 



2.4-1- Strategies I-III: axial/radial focusing Figures|2|l-ni show the spatial variation 
of the departure from harmonicity of Strategies I-III. For Strategy I, II and III the 
departure from harmonicity averaged over a sphere of radius 0.25 a is 0.058, 0.007 
and 0.009 respectively. The radially converging lens of Strategy II, and the axially 
converging lens of Strategy III are almost an order of magnitude more harmonic than 
the single coil radially converging lens of Strategy I. 



2-4-2. Strategies IV-VI: isotropic 3D focusing The isotropic 3D lens of Strategy IV 
rapidly becomes anharmonic as one moves axially away from the lens centre because of 
the axial magnetic field zero illustrated in Figure^IV. In both Strategy IV and V the 
lack of axial symmetry means that there are first-order terms in the field magnitude, 
resulting in a gradient which affects the centre-of-mass motion of the atomic cloud. 
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Figure 3. The departure from harmonicity e(R, Z) of the six different lens 
strategies of Fig. Q is considered for displacements of up to 0.25 coil radii. 
Averaged over a sphere of 0.25 times the coil radius, the departures from 
harmonicity for the six configurations are: 0.058, 0.007, 0.009, 0.385, 1.347 and 
0.148 respectively. 



The main problem with Strategy V is that although both the axial and radial 
curvatures are equal at the lens centre, they vary rapidly with position. The lens- 
centre curvatures are also 8.4 times weaker than the axial curvature at the centre 
of a single coil, leading to longer duration magnetic lensing impulses. The cloud of 
atoms therefore experiences the lens' axial anharmonicities for a greater period, again 
compromising the quality of the focus. 

Figurcs|31lV-VI show the spatial dependence of the departure from harmonicity of 
Strategies IV- VI. For Strategy IV, V and VI the departure from harmonicity averaged 
over a sphere of radius 0.25 a is 1.347, 0.385 and 0.148 respectively. The IP lens of 
Strategy VI is thus significantly more harmonic than the two coil lens of Strategy IV, 
which is in turn considerably better than the axially offset single coil lens of Strategy 
V. 



3. Imaging and minimum cloud size 

The separable equations of motion for a lens that is harmonic in 3D (equations 
©,©) allow us to consider motion in each cartesian dimension as a separate simple 
harmonic equation. It is useful to employ the ABCD-matrix formulation used widely 
in geometrical optics. The position and velocity of an atom along a given Cartesian 
axis, say x, is written as a 2-component vector, and the final and initial vectors are 
related via the equation: 



x { \f A B\f an 
v Xt I \ C V \ v X; 



(16) 
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To simplify the effects of gravity we perform the calculations in a free-falling frame 
of reference. In this frame the free evolution of the cloud is an isotropic expansion, 
described by an ABCD matrix. The influence of a converging or diverging magnetic 
lens can also be described by ABCD matrices, as outlined in |18| : 

M 1 (t)=f; *),M 2 =( COSWT ^ sinuJT ),M 3 =( C ° sh ^ T ^ si f WT V(17) 
\Q 1 J V — uismuiT cos lot J V losiwiut coshwr / 

Matrix M% is the translation matrix for a duration t\ M 2 is the matrix for a converging 
lens of strength to applied for a duration r; M3 is the matrix for a diverging lens of 
strength lo applied for a duration r. It is interesting to note that the sinusoidal 
(exponential) path taken by atoms inside a converging (diverging) magnetic lens is in 
stark contrast to the linear propagation of light rays in an optical lens. 

By multiplying these matrices together, we arrive at the final ABCD system 
matrix. An image (i.e. a one-to-one map of position between the initial and final 
cloud) is formed if the condition B = is maintained. In this case the spatial 
magnification A is the inverse of the velocity magnification D. This spatial compression 
and concomitant velocity spread is a manifestation of Liouville's theorem. The 
theorem states that phase-space density is conserved in a Hamiltonian system. Time- 
dependent Stern-Gerlach forces satisfy the criteria for Liouville's theorem to be valid 
|29| . The cloud extent along x in a given plane is given by: 

a 2 Xf =(Aa Xt ) 2 + (Ba Vx f, (18) 

where o~ Xi is the initial position standard deviation and a Vx . is the initial velocity 
standard deviation. An image is formed for the condition B = 0, but the smallest 
cloud size does not necessarily occur in the same plane. For a single lens system, the 
minimum cloud size occurs very close to the imaging plane. For multi-lens systems, 
the image plane and the minimum cloud size do not necessarily correspond. 



3.1. Imaging solutions 

We note here that the 'thick' converging lens M2 is identical to a thin lens of strength 
C(lo,t) = —lositi(lot) (i.e. the original C entry of the M2 ABCD matrix), pre- and 
post-multiplied by a translation matrix with duration r'/2: 

Ml = ( T 'C')(c 1 ) ( T 'i 2 ) " 

The pulse width r' is defined as 

// n 2 LOT . „. 

t'(lo,t) = -tan—, (20) 

LO 2 

and the notation of primes is used to denote times in the 'thin' lens representation. 
This means that we can use many of the simplicities of 'thin' lens optics, even if we are 
in fact dealing with the more accurate 'thick' lensing behaviour. The effective 'thin 
lens' duration of the pulse r' differs from the actual pulse duration r, but otherwise 
the treatments are identical. In the limit of a short, strong pulse lot — > 0, we find 
that t' — > r. If we wish to consider the diverging lens M3, we merely make the 
transformation lo —> iuj in equation (|19f) - i.e. C = lo sinh(wr) and t' — — tanh 

We can model a single lens system by having a translation of Mi(t[), where 
t[ = t\ + r'/2, followed by a thin lens of strength C, followed by a translation of 
M\ (t'2), where t' 2 — t 2 + t'/2. The physical duration of the focusing is T = t\ + 1 2 + r, 
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however the thin lens system has a total time T" = t\ + t 2 + t' = T — r + r' . For a 
single lens system, the condition B = is met for the system ABCV matrix if we have: 

where 

A = ^, (22) 

yielding a magnification (A — 1) /A. This formalism provides a useful way of designing a 
lens system and investigating its focusing properties. For a converging lens, equation 
becomes: 

WT ' SW= A(1^A)- (23) 
We consider an experimental situation where we fix the total time T = 212 ms and 
the maximum Amp-turns at NI = 10, 000 Amps. The geometry of the lens then fixes 
the maximum strength of the converging lens via u> r 2 = /ib^/to. We can now solve 
to find an analytic result for X(lj,t), which is illustrated in figure 0] for a single coil 
lens of radius 5 cm (uj r — 69.6 rads~ x from equation (|1U|) . The X(lu,t) parameter is 
maximised (and the magnification (1 — A) /A is minimised) when: 

1 -w(T-r) cot wt= 0, (24) 
which has the solution A = 0.929 at r = 23.7 ms. This corresponds to a reduction in 
the atomic cloud size by a factor of —13.1 This is achieved when the pulse duration 
r is from time t = T — r to t = T, i.e. the lens pulse ends at the time of focus. Such 
focusing in three dimensions would increase the cloud density by more than 3 orders 
of magnitude. For a lens placed later in time, the magnetic pulse would not have 
finished at the predicted focal time T, resulting in an increase in cloud size at time T. 




Figure 4. By fixing the strength of a 5 cm Strategy I radially converging lens, 
uj r = 69.6rads — , and total experimental focusing time T = 212 ms, we can vary 
the lens pulse width, t, and find the focusing parameter A(r) via equation 1231 . 
The impulse r has a minimum of 3.9 ms at A = 1/2, and A is symmetric about 
this point. Also shown (dashed) is the result obtained if one makes the strong, 
short pulse approximation lut — > 0, leading to the simplification sinur pa lot in 
equation 1231 - resulting unsurprisingly in a divergence for large pulse durations. 



The above analysis would seem to suggest that the optimum strategy for achieving 
the smallest cloud size would be to construct a lens with a short, strong pulse cut — > 0, 
and use the latest possible pulse time A — > 1. However, experimental constraints and 
lens aberrations alter the above conclusion. 
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4. Investigating various coil configurations for pulsed focusing 

4-1. Methodology 

The ABCT> matrix formalism outlined above is a useful starting point for studying 
pulsed magnetic focusing. However, this formalism ignores magnetic aberrations 
arising due to the departure of the real potential from an ideal parabolic spatial 
dependence (figure |3J). For a non-parabolic potential the change in position and 
velocity which occur during lensing must be calculated numerically. Here, for the first 
time to our knowledge, we test the 'perfect' atomic lens approximations by performing 
numerical focusing simulations. The cloud and its motion are treated classically, 
and for the atomic densities encountered in the expanding cloud, the collision rate is 
negligible. The atoms travel on ballistic trajectories, except when a magnetic impulse 
is applied, in which case the full Stern-Gerlach force is included in the numerical 
integration. For the sake of definiteness we chose to investigate the focusing of 85 Rb 
atoms in an atomic fountain launched vertically through a height of 22 cm (which 
corresponds to a flight-time of 212 ms), a height which is of interest experimentally. 
The atoms come to rest at the apex of their trajectories where they could be used for 
further experiments or loaded into a dipole trap. The effects of gravity were included, 
but these effects on the quality of focus were found to be negligible for the parameters 
used in these simulations. 

The approach adopted was a numerical simulation, in which the trajectories 
of typically 500 atoms are followed. The initial velocity and position probability 
distributions are isotropically Gaussian for each Cartesian direction. The standard 
deviation of position was chosen to be a value typical of experiments at 0.4 mm 0], 
and the velocity distribution corresponding to a typical launch temperature obtained 
with Rb in moving molasses, namely T = 20 /iK. These simulations facilitate the 
calculation of statistically relevant quantities, such as the standard deviation of the 
time-dependent size of the atomic cloud. In all the simulations the maximum current 
value in any coil was limited to 10, 000 A. 

4-2. Strategies 1-111: axial /radial focusing 

To illustrate our methodology we discuss the focusing properties of a single, circular, 
current-carrying coil (Strategy I). Figure[S](a) shows the evolution of the x component 
of a launched cloud of 500 atoms subject to a radially converging lens constructed from 
a single 5 cm radius coil. The impulse is applied half-way in time, and the length of 
the impulse is chosen to reverse the transverse velocity, as can be seen from the change 
in sign of the gradient after t = 106 ms. For this case, the ABCT> matrix predicts a 
radial focus with magnification — 1, which is in excellent agreement with the numerical 
simulation using a parabolic lens (using the strength of equation jlOjl). The vertical 
line at t = 212 ms corresponds to the imaging time. For real coils, it is seen that 
the focussed cloud image is significantly larger than the initial cloud. The aberration 
worsens as the coil radius decreases. Note that although focusing is seen in the radial 
direction, defocusing is seen in the axial direction due to the opposite sign of the 
magnetic field curvature. 

Figure (b) contains an analysis of the cloud in terms of shells of different radii 
measured from the centre of the coil; atoms further from the centre are not focussed 
as tightly, and also focus earlier in time. As the ratio between cloud extent and coil 
radius decreases, the departure of the field from the parabolic approximation becomes 
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Figure 5. (a) Simulation of 500 atoms going through a 10, 000 Amp-turn strategy 
I lens. The solid, dashed and dotted lines correspond to a parabolic approximation 
lens, a 5 cm radius lens, and a 3 cm radius lens respectively. The duration and 
timing of the impulse is chosen using the ABCD matrix formalism, bringing the 
atoms to a focus 22 cm above their launch height, (b) A shell plot of 1000 
simulated atoms passing through a 3 cm radius strategy I lens. The distances 
from the coil centre are 0-10% of the coil radius, 10-20%, through to 50-60%. 
Atoms further from the centre are not focussed as well, and the focusing occurs 
at earlier times. Both of these factors degrade the image quality and size. 



less significant. Therefore one method to reduce the aberrations experienced by the 
atoms is to increase the coil radius, or decrease the atomic cloud temperature. 

We turn our attention to obtaining the minimum cloud size, equation (|18fl . by 
investigating the effect of the A parameter. As discussed in section l3~Tl a thin parabolic 
lens produces the smallest cloud size when A — > 1, i.e. the pulse is applied as late in 
time as possible. This is a manifestation of Liouvillc's theorem - a more compact 
spatial extent can be generated at the expense of a larger velocity spread. 

Figure [S] (a) shows simulations of radial focusing. The radial cloud expansion 
factor, a x /a Xil is plotted as a function of the imaging parameter A. It can be 
seen that the smallest cloud size for a parabolic lens occurs when one waits as 
long as possible before focusing, i.e. A is as close as possible to 1 (limited by the 
solution of equation J22J0- This means that the magnetic impulse ends at the desired 
focusing time. For parabolic fits to the 3 and 5 cm radius single coils (Strategy I) the 
maximum values of A are 0.968 and 0.929, corresponding to negative magnifications, 
M = (1 - A)/A, of 1/30.1 and 1/13.1 respectively. 

Along with the parabolic case, figure (a) shows numerical simulations for 500 
atoms passing through 3 and 5 cm radius coils, for both Strategies I (single coil) and II 
(two coils) . As expected the 5 cm lens better approximates a parabola. Compared to 
the ABCD matrix result there is a marked difference in the behaviour of the minimum 
radially-focussed cloud-size for fields from real coils - the value of A at which the 
ABCD minimum is obtained is dominated by aberrations in the magnetic field. The 
ABCD matrix approach does not provide an adequate description of pulsed magnetic 
focusing when one considers the entire atomic cloud. 

In figure H3 (b) the radial expansion factor of a 5 cm single coil (Strategy I) is 
plotted against time for values of A varying from 0.2 to 0.8. 

The easiest way to reduce aberrations appears to be the use of a very large coil 
radius. Unfortunately the curvature of the field decreases with the cube of the coil 
radius (equation HlOjl ). which necessitates longer pulse durations for larger radius 
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Figure 6. (a) The radial expansion factor, at T = 212 ms, for radially converging 
lenses is plotted against A. The solid line without symbols shows the parabolic 
lens solution, the solid (dotted) lines with symbols shows the result of atoms 
passing through Strategy I (II) lenses. Circular (square) symbols are used for 
the 3 cm (5 cm) radius lenses. The Strategy II lens is more parabolic, allowing a 
smaller radial cloud size to be achieved: a cloud image 0.75 times the original size 
occurs for A = 0.65 using a 5 cm lens, (b) For a 5 cm Strategy I lens, the radial 
expansion factor is plotted against time for A values 0.2, 0.4, 0.6 and 0.8. 



coils. This increase in pulse duration reduces the maximum value of A that can be 
used and therefore also limits the minimum cloud size. The aberrations can only be 
further reduced by increasing the current-turns, something which has experimental 
limitations. 

The aberrations associated with a real coil dramatically affect the strategy for 
achieving a radially compact cloud. The impulse has to be applied significantly earlier 
than an _4,BC2?-matrix analysis would suggest. However, for experimentally-realistic 
parameters it is seen that it is possible to achieve a final radial cloud size that is smaller 
than the initial size. Strategy II (two coils of radius a carrying the same current J, 
separated by S — 0.58), clearly has a better performance than Strategy I (single coil), 
as seen in figure (a). Most importantly, in the 5 cm case, the radial extent of the 
final cloud is 0.75 times the original radial extent of the cloud. 

Strategy III produces an axially converging/radially diverging lens with a high 
level of harmonicity (figure 01 comparable to that of Strategy II. We omit our results 
for axial focusing in this paper, as it will be revisited in the context of alternate- 
gradient focusing in a future publication |19|. 

4-3. Strategies IV- VI: isotropic 3D focusing 

Isotropic 3D focusing can be achieved using two coils with differing currents (Strategy 
IV) . For realistic experimental parameters, the numerical simulations showed that the 
aberrations in the lens smeared out any focusing. However, for unrealistically large 
lens radii and large currents (e.g. 15 cm and 200, 000 Amp-turns) it is possible to 
achieve 3D focusing. 

As discussed in section ^. 3.31 a single coil can be made to have isotropic curvature 
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(Strategy V). At z — ±y2/7a, the axial and radial curvatures are equal, and 
the gradient of the field is non-zero. A numerical simulation was performed for a 
launched cloud, with an impulse applied when the cloud's centre of mass reached a 
distance z = +^2/7a from the centre of a single coil. Due to the large departure 
from harmonicity for the experimentally realistic parameters we used, aberrations 
dominated and focusing was not observed. 

The baseball lens (Strategy VI) yielded the best isotropic 3D lens. Figure 0(a) 
shows the temporal evolution of the volume expansion factor, (a x o~yO- z /o~ Xi o- yi o~ Zi ), 
for a launched cloud subject to a focusing pulse from a baseball lens. Five different 
values of A are depicted, from 0.3 to 0.7 in steps of 0.1. The bias coils have radii 
of a = 5 cm, separation S — 15, and current NI = 3082 A; the baseball has sides 
of length W = L = 2 and carries a 10, 000 A current. The minimum cloud size is 
obtained when A = 0.3, and represents a 31.2 increase in cloud volume at the focal 
time, T = 212 ms. This is to be contrasted with the 13, 000 increase in cloud volume if 
no magnetic lens were used. It is interesting to note that for the parameters we have 
simulated, the results of a 'pure' Ioffe-Pritchard lens (W = 2, L — > oo and 5 = 2) are 
almost identical to the baseball coil. Surprisingly the baseball lens performs better, 
and produces smaller cloud sizes at the focus. 



(a) 




(b) 
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Figure 7. (a) A simulation of 500 atoms sent through an isotropic baseball coil 
lens (Strategy VI). The ratio of the cloud volume to the initial volume is plotted 
as a function of time for values of A ranging from 0.3 to 0.7. (b) A close-up of the 
simulations near the t = T imaging time (vertical line). 



5. Discussion and Conclusion 

We have highlighted the limiting size of a focussed launched cold cloud of weak- 
field seeking atoms for various pulsed magnetic focusing techniques. The ABCT) 
matrix formalism is convenient for giving an estimate as to the parameters needed 
for magnetic focusing, but does not contain the departure of the potential experienced 
by atoms from a perfect parabolic dependence for fields produced by real coils (and 
bars). In this work we have shown how important it is to consider these aberrations 
as they drastically alter the results. We have identified the origin of these aberrations, 
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and described techniques for minimising them. 

We have discussed single coil lenses (Strategy I) , as well as five novel atomic lenses 
strategies (II- VI) , and tested their aberrations both analytically and numerically. For 
all of the Strategies we found that our analytic results for aberrations (section 3) tied 
in well with the numerical simulations of section 4. It was demonstrated in Strategies 
II and III that a 'doublet' radially (axially) focusing lens formed from two coils with 
relative separation S — 0.58 (2.63) provided much tighter focusing than the single-coil 
lens of Strategy I. Amongst the isotropic 3D lenses we found that the baseball lens 
(Strategy VI) was superior to the two coil lens of Strategy IV, which was in turn 
considerably better than the axially offset single coil lens of Strategy V. Of the single- 
impulse lenses, the baseball lens offers the best possibilities for isotropically focusing 
a cloud of weak-field-seeking atoms in 3D. Experiments to test these predictions are 
underway in our laboratory. 

It should be noted that in section 4 we have used the rms radius of a cloud of 
atoms to measure how tightly the entire atomic cloud is focused. By only considering 
a low- velocity fraction of the atomic distribution, even fractions as large as 50%, it is 
possibly to reduce drastically the rms focal spot size, increasing the atomic density by 
orders of magnitude. We will address this complex issue in more detail in our future 
publication dealing with alternate-gradient focusing [T$] . 
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